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ABSTRACT 


A review or the dual-phase power system was made. This 
study focused on the multi-component nozzle of this dual- 
phase system. First, an existing computer code predicting 
the nozzle performance was updated, and second a series of 
experimental tests on a variable area, two-dimensional 


nozzle was performed to verify the computer code. 
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The dual-phase nozzle is a key element in the dual-phase 
engine concept. This nozzle will be studied in detail, but 
first a review of the dual-phase cycle will be carried out 
from information obtained from References 1 and 2. 

Reference 1 describes the dual-phase engine. It isa 
new concept which operates on a mixture of two fluids or two 
phases of one fluid. This cycle employs a high-torque/low- 
rpm impulse turbine which eliminates the requirement for 
speed reduction needed in a conventional steam turbine. 

Maas allows for direct drive of a ship's propeller with no 
gear reduction. This is ideally suited for marine propul- 
Sion, reducing weight, noise, and volume. 

The dual-phase concept is a modified Rankine cycle with 
subtle differences of extreme importance. In a normal steam 
turbine the working fluid enters a turbine through a nozzle 
where the kinetic energy of the working fluid is converted 


to a mechanical form. The dual-phase system introduces a 


tt 
is 


meeend fluid prior to entry into the nozzle. This fluid is of 


{ 


higher vapor pressure than the steam and therefore remains 
M@@eene liquid state throughout the cycle. Section A will 
describe this two-component cycle while section B will do 


the same for a single-component system. The dual-phase 





system can be divided into the two groups illustrated in 
Figure 1. Single-component flow can be further divided into 


three categories. 


A. TWO-COMPONENT 

A two-component mixture is one in which the low vapor 
pressure liquid and a high vapor pressure liquid are of dif- 
ferent chemical compounds. Some fluid combinations which 
have been considered are steam-krytox, steam-caloria, steam- 
lead, bismuth eutectic, and dow-therminol. A schematic flow 
diagram and process representation on the T-S diagram are 
shown in Figures 2 and 3 for the two-phase engine cycle using 
a "two-component" mixture. The liquid phase is lithium car- 
bonate and vapor phase is steam. To illustrate the overall 
advantages of the two-phase engine cycle a discussion on the 
theory of operation will be presented using a two-component 
mixture and Figures 2 and 3. 

The major component of the dual-phase system is the nozzle. 
Pemixer area is located prior to the nozzle inlet. A high 
meper pressure liquid is placed in contact with the low vapor 
pressure liquid in this area. A high pressure vapor liquid 
mixture is formed. Since the temperature of the liquid is 
greater than temverature of the water, heat is transferred 
tO the water causing it to vaporize to point 1. Figure 4 
lllustrates the temperature and state point from the inlet 


t= 


mene Mixture area to the nozzle exit. This mixture is 
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expanded to a low pressure at the nozzle exit. The expand- 
ing vapor transfers momentum to the liquid droplets, while 
the droplets transfer heat to the vapor. This tends to ap- 
aproach isothermal expansion through the nozzle, point 2. 
After leaving the nozzle, the liquid droplets are separated 
from the vapor by a rotary separator located on the turbine. 
The rotating drum turns at close to the nozzle's exit vel- 
ocity. Hence there is very little friction loss between the 
separated liquid and the drum wall. The centrifugal accel- 
eration produces a very clean separation of vapor and liquid. 
The liquid traveling with the rim of the rotating drum trans- 
fers axlally through holes in the separator dise and subse- 
@uentily enters the liquid turbine proper. The kinetic 


energy of the liquid is converted to shaft horsepower by the 


Hrqguid turbine. Figure 5 illustrates the dual-phase liquid 
impulse turbine assembly. The vapor remains superheated as 
pmeonsequence of the heat transfer from the liquid. The 


vapor flows from the separator, point 3, through a regenera- 
tor where heat is transferred to the condensate. The vapor 
is condensed, point 5; pumped to nozzle pressure, point 6; 
and passed through the regenerator for heating. Heat is 
added between point 6 and point 1 in two methods. The re- 
generator adds heat to the condensate by using the steam 
trom the turbine rotor, point 7. The remainder of heat, 
point 7 to 1, is added by the heated liquid mixed with the 
Condensate in the mixer. The water is vaporized by direct- 


@emeact heat transfer. 
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Ty 
ar 
MIXER 
Ty 
where T,. = Temp of liquid entering mixer 
a 
T 4 = Temp of water entering mixer 
Ih 
ai i 
Ts = Temp of steam leaving the nozzle 
7 = Temp of liquid leaving the nozzle 
ay. 
Figure &. Temperature & State Point Diagram 
MeoGecnenliixsure ¢ Nozzle 
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immeehewmNeZ2"e Most Of the thermal energy of the steam 
if converted to kinetic energy of the liquid droplets. 
This acceleration of the liquid by the vapor in the two- 
phase nozzle provides the kinetic energy to drive the 
Jiguid impulse turbine. The liquid velocities involved 
are relatively low as compared to velocity of the vapor. 
Meme, the output of the impulse hydraulic turbine will be 
aon torque/low rpm. This conversion of the liquid kinetic 
energy to shaft power at high torque with low rpm appears 


memiave direct application to naval propulsion. 


bee ONE-COMPONENT 

A one-component system is one in which the working fluid 
1s of the same chemical compound. One of the simplest dual- 
phase one-component systems is illustrated in Figures 6 and 
7. The working fluid is heated to saturation temperature by 
some type of heat source. Heat sources applicable to this 
case are geothermal power plants, engine exhaust, industrial 
waste-heat recovery, and bottoming cycles for steam and gas 
Memaine Plants. This working fluid, at saturated liquid con- 
ditions, with small amounts of vapor is expanded through a 
two-pnase nozzle. As the expansion process takes place ,the 
liquid partially vaporizes and accelerates the remaining 
liquid phase in the nozzle. The dual-phase mixture enters 


the rotary separator and the same process occurs as mentioned 


ae 


in section A. Since the liquid phase is of a much higher 
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density than the vapor, the liquid velocities are relatively 
low aS compared to the velocity of the vapor. Thus the out- 
Peeeor the turbine will again have high torque at low rpm. 


After leaving the turbine the vapor mixture 1S condensed and 


the condensate is pumped back to the heat source. The cycle 
is shown on a T-S diagram in Figure 6. The state points are 
numbered to correspond with Figure 5. The dual-phase nozzle 


Sepenclon takes the fluid from a saturated liquid, point 1, 
to a dual-phase flow, point 2. The flow is decelerated in 
the rotor; condensed, point 4; and pumped back to nozzle 
maret pressure at point 5. The liquid is then reheated by 
ieee source fluid to point l. 

Another application of the one-component two-phase cycle 
HM the wet-to-dry cycle. If the initial temperature of the 
@em<ine fluid is sufficiently high and the saturation curve 
has a positive saturated liquid slope the working fluid can 
Be expanded to dry vapor. Figure 8 is the T-S diagram for 
a wet-to-dry cycle. The fluid 1S expanded from saturated 
liquid at point i to saturated vapor at point 2. The vapor 
drives animpulse rotor and leaves the rotor slightly super- 


Meeted at point 3. The vapor is condensed to point 4 and 


Pumped back to the nozzle inlet pressure at point 5. 


C. ADVANTAGES 
The advantage of the dual-phase cycle with respect to 


Meeeine appiication is the ability to achieve low shaft speed 


29 





TEMPERATURE, °C 


300 


250 


200 


150 


100 


50 


A SATURATION 







CURVE 
SOURCE 
FLUID 
NOZZLE 
HEATER 
CON DENSER : 
bs) | ROTOR 
PUMP ; 
ENTROPY 
Figure 8 Wet-to-Dry T-S Diagram 





in a small turbine engine. For example a steam turbine of 
150-kw shaft horsepower using the temperature given in 
Figure 2 would have a speed of about 60,000 rpm. A compar- 
able two-phase turbine would operate at approximately 10,000 
rpm. There is also an efficiency advantage. At temperature 
corresponding to Figure 2, a steam Rankine cycle would have 
an efficiency of approximately 28% where-as the two-phase 
eycle would have an efficiency of 37%. This is assuming equal 
turbine efficiencies. The two-phase cycle also allows for 
control of turbine speed because the vapor/liquid mixture 
ratio can be varied to change the nozzle exit velocity. This 
is a capability unavailable in a conventional steam turbine. 

Both of the dual-phase concepts can be thought of as a 
form of a regenerated Rankine cycle. The dual-phase cycle 
by control of liquid/vapor mixture ratio enhances the overall 
power system controlability. The T-S relationship for a dual- 
phase two-component engine cycle compared to a Rankine cycle 
Memshown in Figure 9. 

Two design studies References land 2, have shown potential 
advantages in the two-phase engine cycle as compared to the 
Semventionel Rankine cycle for marine propulsion. The fol-~ 


lowing advantages were noted: 


[v4 


1. High Efficiency - Full load output power verformance 
gains ranging from 20 to 50 percent was found. 


2. Direct Drive - Direct drive at speeds ranging from 


J0-4500 rpm was found possible with a single-stage turbine. 
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3. Reduced Volume - Volume reduction of 30 percent 
were estimated. 

4, High Part-Load Efficiency - Variable mass ratio 
enabled part-load (cruise) efficiency gain of as much as 


myo percent. 
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Peover s ee NOZZLE THEORY 


The flow phenomenon of a two-phase mixture has been ana- 
lyzed in Reference 1. It is repeated as follows. The problem 
is illustrated in Figure 10. A spatially uniform two-component 
mixture of liquid drops and gas enters a nozzle at high pres- 
sure and low velocity and expands to low pressure and high 
velocity. The objective of the analysis is to determine, for 
a specified pressure the drop diameter D and the tempera- 
tures Ds and T,, velocities Us and V,; and flow rates m, and 


So 


my of the gas and liquid phases, respectively, at each station 
Mmeene nozzle given the initial values of D, T_, T;, V_,V,, 
a re 


the total flow rate, and the properties of the fluids. 
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Figure 10. Dual-Phase Nozzle Flow 
Geometry & Nomenclature 
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The five relations employed to compute the five unknowns 
Be ve Th > Vs and V,, are (1) the momentum equation for the 
mixture, (2) the energy equation for the mixture, (3) the 
drop drag equation, (4) the drop heat transfer equation, and 
Moeeethe drop breakup criterion. Solubility and vapor pres- 


sure relations provide the flow rate ratio m/m, . 
oS 


Swe SOoUMPTIONS 
The assumptions employed in the two-component analysis 
eee as follows: 
feline liquid is wniftormly dispersed as spherical drops 
all of the same diameter. 
gs) Ihe drops break up to limit the Weber number to 6. 
fee ihere are no external forces acting on the two-phase 
mixture other than pressure and wall shear, and 
there is no heat transfer to or from the mixture. 
4, The flow is one-dimensional. 
mer lne Grops are large enough for the surface curvature 
to have negligible effect on the vapor pressure of 


the liquid and for the surface energy to be negligible. 


6. The drops are isothermal. 
7. The gas mixture obeys the additive-pressure law. 
8. The partial pressure of the predominantly liquid 


component in the liquid is given by Raoult's Law. 


tO 


The concentration of the predominantly gaseous com- 


ponent in the liquid is given by Henry's Law. 
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10. The volume of the liquid solution is equal to the 
sum of the volumes of the pure liguids. 
Assumption 1 restricts the analysis to nozzles having 
Seacially uniform injection of the liquid into the gas and 
operating at gas-to-liquid volume ratios greater than unity. 


Assumption 2, the drop breakup criterion, states that drop 


diameter is limited to a value D for which We = Ve D/2 
wo. nus 
Doe eee 
piVs (1) 
Peere 9. is the gas density, We fe ene slips velocity Ve - Vis 
© 
and Domi surprece tencaon. Ihe form of Eq. (1) is 


physically reasonable in that the Weber number is propor- 
tional to the ratio of stagnation pressure 0 Ve Be to sur- 
Face tension pressure ¥V/D. Hence, a drop would be expected 
emeerceten and breakup at a sufficiently high value of W,. 
This has been verified experimentally and the critical Weber 
Mmemoer found to be 6, within a factor of about two. An 
@@atcional restriction is that for actual breakup to occur, 
the time spent at a Weber number exceeding 6 must be longer 
Seemethne natural period of oscillation of the drop, 

mp, D°/ ite where OF momene Gensicy Of the liquid. 

This requirement is met only in two-phase nozzles longer 
than about 10 in. and Assumption 2 may cause the analysis to 


Peeeeeotimate the exit velocity by increasing amounts as the 


nozZZle length decreases below 10 in. 


Z7 





Assumption 3 excludes magnetohydrodynamic and mechanical 
body forces. The exclusion of wall heat transfer is correct 
for the insulated nozzles of interest for power systems. [In 
addition, the relatively high velocity results in short 
residence lines in the nozzle proper. 

Assumption 4 is closely met in practical nozzles since 
good performance requires small wall angles, large throat 
Mmernus Of curvature, and uniformly distributed injection of 
the fluids at the nozzle entrance. 


Assumption 5 is valid for the drop sizes of 0.001 to 


tJ 


MoO in. produced by the Eq. (1) breakup criterion. Assump- 
fem) o is valid because of the rapid internal circulation 
mmearops. Assumption 7 introduces negligible error in most 
cases of practical interest since the vapor pressure of the 
liquid is small and needs only to be evaluated approximately. 
Moesumptions 8, 9, and 10 are either valid, or cause little 


peeeee. tOr Fluids of low miscibility, wnich are the fluids of 


imeerest. 


See Ri VATION OF EQUATIONS FOR FREE-STREAM FLOW 
ime Continuity 


Referring to Figure 1, the nozzle flow area A is 


oO 
a) 


equal to the gas flow area Ht ae bius the liquid ftiow 


© 


, ee 
area m, /o,V,- Thus 





i = ee 
A= mm (Fy, a7) 


where ris the mass mixture ratio m,/m_. 
— S 
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2. Momentum 
By Assumption 3, the only force acting on the free- 
Stream flow is that due to the pressure gradient. If M is 
the momentum flux at flow area A, the change in momentum 


flux across pressure increment dp is 


dM = —Adp (3) 


Maes Memenrumeslix Came pe Written as the sum of the 


Mementum fluxes of the gas and liquid. Thus, 


Af = mV, + mV; ee 


If the flow were allowed to continue at constant 
pressure, Vg su gyal is would become equal to each other at the 
mass-weighted mean velocity V. Since for this process, 


@ee- 0, the value of V is given by 


(m, + rit:)V = mV,+ mV; C3) 
ON Gy 
" Vi +rV 
va aa oe 


Thus, the momentum flux can be written 


M =(m, + m)V ie 


pence m_ + ny is constant, the change in momentum flux is 
g | 


; ; eee (3) 
dM = (m, + m;)dV 


meee tuting Eqs. (8) and (2) into Eq. (3), dV can be written 
a 1 1 a al d 
dv = eZIWAS prV, 7 (9) 


to 





The slip ratio is defined as 
2 = V/V = (V,-Vi)/V eee 


This equation can be combined with Eq. (6) to give Ys and 


fein terms of V: 








ie 
8 No Ls 
ve= (1+ 755) =e Crip 
Ss 7 
v= (1 r7)¥= Gi) 


The gas density can be expressed as 


3 
ps = W,p/RT, ed 


where we is the effective molecular weight of the gas mixture 
ies 18 the universal gas constant. Eq. (13) is the defini- 
tion of the effective molecular weight ee which is the 
quantity that gives the actual gas density when substituted 
momeeqg. (13). 

eerie ees. (tl) /= (13) into Eq. (9), the dif- 


ferential momentum equation is 


— = — Z. RT, Mik 
BO ave a7 (ans a on ay (14) 


The quantities a and b are slowly varying because s is 
Beeecalily only 0.1 to 0.3 and slowly varying. The quanti- 


ties r, ee Ne and 9, are also slowly varying. Integrating 





Pam (l4) over a pressure increment Ap, for which a, b, r, 


Eee We and p,; are constant to within the desired accuracy, 


the change in V’ is 


cs Z RT,,, p+Ap/2 . fmAp 
AV = 1 + fe ae log. p— Ap/2 — a 


C15) 


All quantities other than pressure can be taken outside the 


integral and evaluated at their mean values (denoted by sub 





Beripc m) corresponding to the mid-interval pressure p. Thus, 
= 2 
ee ee 
ol ae 
RT pe 3 d pee 
2 
x be ee eesenaion he (16) 
onW’s,, Ap P OmP te Ap 
re as a a 
(16) 


os ae: RT r 
a Pace, (ena eee 
ey Se = bor) CET 


Equation (17) is the final form of the momentum equation. 
a -£ 


The enthalpy change of the mixture between state 1 


(the beginning of pressure interval Ap) and state 2 (the 


(D 


nd 





of the interval) can be evaluated in two steps: (1) phase 


change at (0; a al and (2) change to Pn» Bone Th? at 


ae Ly 
fixed composition. 
The enthalpy change for step 1 is 


enthalpy required 


AH eaMounic 1O fe to vaporize and 
Ee vaporized Xx heat unit mass of 

ee eomra, to ) 
L el 
amount of B enthalpy required 

rr vaporized x to vaporize and 


Mest mass “OL 


Boom eeom qT wel ab 
zt 8) 


kinetic energy 


amount of A required to 
+ and B xX accelerate unit 
mass from Ve to 
V 1 
8] 


ey 


4H, = (™Ma,, = Ma, |) [Lc, 3 Cay, (T,, — r,)] 
ze (Ms, = m»,.) [L., a Co, (T,, a Urea See 
an (m,, a2 m,,)(Vi— NOE )/2 


Mepees |, and c are latent heat and specific heat. respectively. 
15 | 


Introducing more compact notation. 


AH, = Am,, (Le, + Ca, 5:T) 


° m.5.V2 
+ Am), (Le, _ Co, 3:7) _ — 


Piewemmenadlpy Change for step 2 1s evaluated from the 


temperature, pressure, and velocity changes, with properties 


Peeemeuated at mean T and o0 for the interval. 


SZ 





AH, 


1 : 
Ty, [¢s., (T,, . T,,) a3 “o (V os V3.)] 


25 iy ey = Wee — 


C20), 
1 oa: 2 
+ 9 ( ts Vi,)] 
AV? 
= Me, (< ro ME gs =") 
Mo Oe 
+n, (61,7: + 22 + =F") (21) 


By Assumption 3, no work is done by the free-stream flow 


and no heat is transferred to it. Hence, 


AH, + SH, =0 (22) 


Mp orhiturine bac. (19) and (21) into Eq. (22) and 


eeeving for aa gives the energy equation for the mixture: 








AV: 
Ars a2! = (6... 47 +22 4 svi) 


Ce. 9 


coe “a 





am,5,V? Am, 
~ Ce > 9 : 


Arty, 
+ my, (Ls, + Co,. 7) 


Cle 


eee Drag 


Although no force other than pressure acts on the 


free-stream flow as a whole, a drag force exists between 





the phases. Hence, a second momentum equation must be 
witten using as the control volume the boundary between 
the phases. 

The two forces acting on each liguid drop are the 
buoyancy due to the pressure gradient and the drag due to 
the relative gas velocity. The sum of these is equal to 
the mass times the acceleration of the drop. Thus, for a 


Single drop. 





dynamic pressure of drag Frontal area 
relative gas flow x @Gectiicient: <x Or dean 
7 volume pressure 2 mass y acceleration 
oe drop gradient drop drop oi vdrop 
or 
] xD? *D* dp 
5 Pr|V.|¥s) Co tae re 
(: oo (24) 
. dV 
= mae vi) 


The absolute value Sign in the first term makes the 


drag force positive when V_ > ie and negative when V_ < V 


g g ne 
Solving Eq. (24) for dv, 
ae Sp,|s|sV?Codx _ dp 
ai 4p,V,D piv: 
C25) 





Poerterentiating Eq. (12), dV, can also be expressed in 


terms of s, r, and V. Thus, 


HS sdr  __ds 
av, = bd + 0 oF eat (26) 


perving for ds, 


@ Seed av 7 sdr__ (1 + +) dVr 


V l+r V e275 


upstate mne sv mpemetg., (25), noting that dV = 


L 
Ba / 2V. using £q. (12), and writing for a finite increment, 


results in, 


b.. (1 + tm) AV? ae (1 + 7,,.) Ap no SmAf 


As = ove bap: V2, 1+ Ton 


(23) 
3p,,, |Sa|SmCo,, (1 + tm) Ax 


4 bapi,D 


This is the drag equation employed when x is specified 
Meme function of o. 

Solving Eq. (28) for A x yields the required alter- 
Mative equation: 


4D b*, pi, AV? 
= a aa. ae jap + ao 
SPs. | Sun | St ye | a= 


+ bn Pt, ve SmAP 
oe Tm 1+ lm = 


Sie 





5. Heat Transfer 

Although no heat is transferred to the mixture as a 
whole, heat transfer exists between the phases. Hence, a 
second energy equation must be written using as the control 
volume the boundary between the phases. 

The work dW done onthe liquid is that due to drag 
by the gas. (Only work done by shear or shaft forces is 
included in dW when writing the First Law for a control 
volume). Multiplying Eq. (24) by the number flow rate of 
drops N = bit, /TD°p, the drag force a Onekhar, Guantity of 
Peouid is 


N m, d : av 
Senator) || scion) pis OY Tae (30) 


mmeework done on the liquid is 


Z 26 do 2) 
—dW = Feds = su (F ae: ecue 


The heat dQ transferred from the liquid is made up of 
two parts: (1) the convective cooling due to the temperature 
difference between the liquid and gas and (2) the evaporative 
cooling due to the latent heat supplied to the liquid vaporized. 


The convective cooling is 


ao: = hAWN(T; — T,)a? 





where h is the heat-transfer coefficient, A 4 = ns 1s the 
surface area of a drop, and dt = dx/V_ is the time required 


to traverse dx. Thus, 


=—dO-z _ 6hmi(T, a T )dz 


Dp.V, C33) 
The evaporative cooling is 
—dQ, = L.dm,, + Lidmy, (34) 


The change is enthalpy of the liquid over the pressure 


inerement dp 1s 


a dp ot) 
aH = ra, ( edt, + 2 + (35) 


Peseta eine Ease ol). (33), (C38), and €35) into the 


Steady-flow energy equation dQ - dW = dH, the result is 


6hm ,sT dz 


Dai: = L.dm,, - L.dm,, =m,c,dT; 


Colles, 


Where 6T = T_-T,. 
g L, 
Pemrene for a fimite anterval, the final form of fhe 


meep heat-transfer equation is 








‘AT, =-—— yaa mi Mi 


1 fon bine) Te aul ma A | 
Ci 


a7 





power meme mmm ge (23)5°(23), and (37) are the five 
Equations that must be solved simultaneously to obtain the 
values of the five dependent variables D, Des Th > Vig and 
ua eset lOonrOrnitne Independent variable p. To carry out 
the solution all quantities in the equations must be expressed 


moecerms of these six variables. 


C. WALL SHEAR AND BOUNDARY LAYER 
For a two-phase nozzle, the momentum flux of the friction- 


less nozzle flow is that given by 


W=mV 


The mean mixture density corresponding to the mean velocity 


Vis 


3 
> 





| 


D 
se 
< 
— 
a 
a 


a8) 


where ee 1s the ratio of gas flow area to liquid flow area 


ep *°.',- 


From the definition of the momentum thickness, the value 
Memo at a4 Station where the nozzle wall radius is She is 
given by 

M — M, = 2ny.8 a = 2ny,p’V26 
(39) 


ferme M- is the momentum flux of the real flow with friction. 


i 
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The skin-friction coefficient can be defined using the same 


quantities as single-phase flow 


Ore 
Or = Taye (40) 
where sm moneyed. Liteweell, be shown that a valid 


C- value can be peovVEded . 
The boundary-layer momentum equation then becomes 


dj = a dx 





I + (8*/8) = x = | 
a| . Be pan.) R, aR | 


where 6* is the displacement thickness, i.e., the distance 
the wall must be moved outward to give the same flow rate 


Memeotten frictionless flow. 


Assuming a $ power velocity profile and no density vari- 


fmeeom, the shape factor 6*/&S is obtained from 


i “EL = (y/3)""] dy 


d | ° (y/sy" [1 — (y/8)'") dy 


5° 





~i| co 


where 56 is the velocity thickness of the boundary layer. 
Noting that dV can be written a¥*/2v, inc chvcmereadt a. V)\= 


d(m,/A), the fimnite-ditference form of Eq. (41) is 





Cn ey ee 
a _ = a ™ (aH = Net ° Yo hs 


nm 


‘ee) 
oO 





Wall shear in homogeneous two-phase flow has been found 
to be equal to that which would exist with pure liquid at 
equal velocity and boundary-layer thickness, multiplied by 


lmmenwetted wall fraction: 





5 (42) 
5 D0 =a 


where aa 1s the skin friction coefficient for liquid at a 


Reynolds number of 


Pepeonvenient relation for Ce as a funet.on OT Ro Ss, the 
i. 
Sshultz-Grunow relation which can be written 


Ee 
C;, ase (logic Rs a A252 


eee "8 


Comparison of Eqs. (40) and (42) shows that C- can be 


written 


Thus, the final form of the boundary-layer momentum 


Eeuation is 





4 Q 





Let V, be the mean velocity of the flow including the 


Bemndary layer. Then, from Eq. (39), 


mn.Vs = M,= mV — 2ny,p'V26 


Hence, employing Eq. (38), the mean exit velocity includ- 


ing the boundary layer is 


= _ oY 
¥o=0(1- =") 


By the definition of the displacement thickness, the flow 


rate is reduced by the throat displacement thickness 46% to 


+ 
( =) 
e _ m l a pA Pee 
ms t UD, 


‘er 


Dee NOZZLE THEORY SUMMARY 

The preceding equations form the basis for the mathematical 
model which is used to predict, based on inlet conditions, the 
exit velocity, and temperature of the mixture. These equations 


See, form the basis for the model which provides the optimum 


nozzle shape given a set of inlet conditions. Some additional 
relationships are, however, required. These are: 
1. Phase properties - to establish the mass ratio, mass 


flow rate ratio of gas to liquid, and the thermal 


conductivity of the mixture. 





woeoiauid drop drags coefficients. 
3. the liquid drop heat transfer coefficients. 
4, Boundary layer momentum thickness and displacement 
thickness . 
fee okin friction coefficient. 
These five additional relationships are developed in detail 


mmereference [1]. 
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ITI. COMPUTER PROGRAM DUAL-PHASE NOZZLE 


The computer program employed in this study is based on 
mmmorram developed by Dr. G. Elliott of the Jet Propulsion 
Laboratory in Pasadena, California. The program was updated 
and converted for use on the Naval Postgraduate School 
computer. The Dual-Phase Two-Component program employs the 
theory in Section II. The program is written in Fortran 
computer language and can be compiled using a Watfiv or For- 
jean LV compiler. 

This program has been utilized in dual-phase nozzle 
analysis and to provide values for comparison with the ex- 
perimental results. To use the computer program the inlet 
GonGcitions have to be specified. The flow conditions are: 
inlet pressure; mass ratio; inlet temperature of the gas and 
giguid; inlet velocity of the gas and liquid; total mass flow 
rate; and nozzle exit pressure. Section III B, shows speci- 
m@easdetails for data input. 

There are two options that can be chosen. The first is 
prescribed pressure-versus-distance option MOP=0. The vres- 
sure profile P(X) is selected corresponding to the adopted 
nozzle contour. If the pressure-versus-distance is used a 
P(X) input table is required. This profile is developed 


from the actual measured pressure values in the experiments. 


U3 





(See Appendix A for sample program.) The second option con- 
icfcnoL dl Optimum nozzle contour option MOP=1. This option 
is useful only when the liquid drop diameter is constant. 

The dual-phase two-component computer program is a struc- 
tured program with thirteen subroutines controlled by a main 
program. This arrangement improved the programming process 
through better organization and programming notation. 

The control point of the dual-phase two-component com- 
puter program is the "main section." It controls the flow 
path and operation of all input data, property tables, and 
calculations. It accomplishes this by calling the thirteen 
Subroutines at the appropriate times, saving wanted data 
in files, and printing out desired information. 

One of the most important subroutines which inputs in- 
BemmatiOn is the "“INTRP™ subroutines. INTRP controls the 
property table inputs. It reads in four two-dimensional 
tables and fourteen one-dimensional tables. These inputs 
Meemene properties of the gas and liquid phases of both 
components of the flow in the nozzle. The subroutine writes 
the values of these tables into a file and retrieves appro- 
Betacte values from that file. INTRP can also interpolate 
for values used throughout the entire program. 

mmeut Of case data 1s controlled by subroutine "Sect 1." 
identification information and case heading information is 
read and printed for each instance. Sect 1 also places the 


Seeeooure vs. distance profile, if specified, in an array. 
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"Sect 2" through "Sect 6" are the subroutines which cal- 
culates the flow data. Sect 2 sets the initial conditions 
indicated in the input, and begins the iterations. Sect 3 
Pemputes initial flow rates of both components; the initial 
area of the nozzle; slip velocity; mean free stream velocity; 
and slip friction. Sect 4 computes the changes in flow par- 
ameters and new distances. It then begins to calculate new 
conditions such as flow rate, temperature, velocities, sur- 
face tension and mean area. Sect 5 is the binary cut con- 
vergence routine and computes mean boundary-layer parameters. 

If a problem is diagnosed in any subroutines and "diagno" 
meecalied, it will print all output parameters calculated to 
that point. It also does the same if there is a convergence 
problem. 

There are two subroutines that output calculated data, 
S@eroutine "Write" and subroutine "Output." Subroutine 
eee” will send output information to the printer for a 
Mamaepaper copy. Subroutine "Output" reads and stores the 
Sa@eput on a file. 

The two-phase two-component program has been written 
with comment statements in the text of the program. These 
will allow for a more understandable and, therefore, a more 
Easily modified program. For specific details on the content 
of these subroutines, see Appendix M. 

The dual-phase nozzle program was tested for correct 


Output. Sample data and results were obtained from 





Sree David Elliott, were inputted into the program, and 
executed. The output was compared to the sample data. The 
program produced duplicate results. 

The program begins by storing fluid property tables and 
reading in all input data. All nozzle inlet conditions are 
Semputed. The program then proceedshalf a pressure step at 
atime. At the middle of each pressure interval, the changes 
in quantities across the interval are computed using the pro- 
perties interpolated from the table for that pressure, and 
for the existing temperature. The change in slip is found 
aft the pressure profile P(X) is specified. At the end of 
each pressure step, the flow conditions are updated and 
initial conditions are determined for the next step. The 
dropsize is reduced at the point when the Weber number ex- 
ceeds six. The flow conditions are printed if the pressure 
m@eamone selected for output. The computation continues until 
the last pressure step has been completed and flow conditions 
at the smallest flow area encountered are printed as the 


moreeat conditions. 


See ROPERTY TABLES 

i weat Capacity of component "A" vapor in BTU/LBM-R 
is a function of temperature and pressure. The two-dimen- 
Sional tables are entered row-wise. At least two cards are 
mecessary to specify a row and at least two rows must be 


entered. 
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@anceiiee (tommat GE12.6) 
eols. 
1-12 temperature (R) 
ie=-2teIe Ott this is the last temp for 
this table 
@a~ra 2?  (tonmat bE12.6) 
cols. 
1-12 pressure (psi) 
13-24 heat capacity (BTU/LBM-R) 
25-36 pressure (psi) 
37-48 heat capacity (BTU/LBM-R) 
49-60 pressure (psi) 
61-72 heat capacity (BTU/LBM-R) 

The maximum entries of temperature are 35 values. for 
each value of temperature, the maximum number of entries of 
pressure and heat capacities are 35 values. Each row of this 
table will be terminated with the pressure and heat capacity 
Saat to lige These two values are not counted in the maxi- 
mum of 35 entries/rows allowed. 

The program shown in Appendix B can be used to determine 
gemma ues of heat capacities. The program structures its 
Output in the format needed for the table input. It uses 
input data obtained from Reference 3. The input data must 
De placed in a two-dimensional table. This table is used 
in the program to interpolate the values needed for output. 


Input data must be formatted as follows: 


Lu] 





Card 1-16 (format 1OF7.4) 
cose 
1-7 temperature (R) 
at this temp the following is entered: 
8-14 heat capacity at .01 P 
Te-21 Heat capacity at .4 Pa 
22-28 heat capacity at .7 Pa 
Zoo tess Capacity vat 1.0 Pa 
34-42 heat capacity at 4.0 Pa 
43-49 heat capacity at 7.0 Pa 
SUS coenecaieeapacity au 10.0 Pa 
97-63 heat capacity at 40.0 Pa 
64-70 heat capacity at 70.0 Pa 
The temperature must be entered with 
increasing value. 
Peeeneat Capacity of component "B" gas, BIU/LBM-R is a 
function of temperature and pressure. This two-dimensional 


table has the same format as part Al above. 


3. Molecular weight of component "A" vapor is a func- 
tion of temperature and pressure. The two-dimensional 
tables are entered row-wise. At least two cards are neces- 


sary to specify a row and at least two rows must be entered. 
Page Jee (tormat 6E12.6) 
‘elon yeas 


1-12 temperature (R) 
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13=24=1.0 if this 1s the last temp for 
this table 
Card 7: (formateonl2eo) 
Cote. 
1-12 pressure (psi) 
13-24 molecular weight 
25-36 pressure (psi) 
37-48 molecular weight 
49-60 pressure (psi) 
61-72 molecular weight 

The maximum number entries of temperature are 35 values. 
For each value of temperature, the maximum entries of pres- 
sure and molecular weight are 35 values. Each row of this 
table wiil be terminated with the pressure and molecular 
weight equal to male These two values are not counted in 
the maximum of 35 entries/rows allowed. 

The program shown in Appendix C can be used to deter- 
mine the values of molecular weight. The program formats 
1tS Output in the format needed for the table input of 
the two-component two-phase computer program. It uses 
input data obtained from Reference 3. The input data 
must be placed in a two-dimensional table. This table 
1s used in the program to interpolate the values needed 
for output. The program used in Figure 8 can only be used 


With ideal gases. Input data must be formatted as follows: 





Card 1-16 (format 10F7.4) 
eos. 
1-7 temperature (R) 
at this temp the following is entered: 
8-14 density at .01 P 
15-21 density at .4 Pa 
22-28 density at .7 Pa 
22-6 Scenics me rae 1.0 Pa 
34-42 density at 4.0 Pa 
43-49 density at 7.0 Pa 
50-56 density at 10.1 Pa 
SJ=ed density at 40:0" Pa 
64-70 density at 70.0 Pa 
The temperatures must be entered with 
increasing value. 

4. Molecular weight of component "B" gas is a function 
of temperature and pressure. This two-dimensional table 
nas the same format as part A3 above. 

9. There are fourteen one-dimensional tables. The 
One-dimensional tables are entered in the following format 
Meae Z(T)): 

Sard Vie Crormat 6E 12.6) 
Gous . 
ee. VE 2 OR 


13-24 Ti-2 


20 





Zoo la = Z 


lon 
am 
J A 
wn 
S 


37-43 Zi-O1 
49-60 Ti 


bl—7 27a and I-11 < TF neue an eh 


Each table is terminated by two consecutive entries of 


5 


Mime 1 .e., Lk 


3 


ioe ae > meeht adjusted in the field). 


The fourteen one-dimensional tables are (in order requested): 


i 


~J 


CAL(T) heat capacity of component A liquid, 
pease OR 

enn) heat capacity of component B liquid, 
Srey EBM OR 

LACT) latent heat of vaporization for com- 
ponent A, BTU/LBM 

Ee cr) latent heat of vaporization for com- 
ponent B, BIU/LBM 

PEO Ll) vapor pressure of component 8 
ROAL(T) density of liquid component A, 

ei) es 

LOL meerdenstey of Iigquid component B, 
INS VBU a dics: 

KAG(T) Eherial cenduectivity of component 
Peecaoeee Uy El HR OR 

KBG(T) thermal conductivity of component 


a idommen o,f l HReOR 


sya 





De linet ye viscosity Of liguid component A, 
LBM/FT HR 

tie some wevicesgcity Of liquid component 8, 
Joy 12 1h Sistas 

Ieee VEAG@oeeviscOsity Of gas component A, 
JL eyOU Aaah jalsst 

13. VIBG(T) viscosity of gas component B, 
Bet ay nik 

ie CLGCh) surface tension of liquid component 
B, DYNE/CM 


Appendix D gives sample property table for input. 


See CASE INPUT 
A blank card must separate the property table from the 
data set decks following. 
Card 1: (format 4A4, A2, 3A4, 7A4, A32, 216) 
eos 
1-18 date 
19-30 case number, may be any 
alphanumeric data 
Si-O 0Mdenrir Leataon 
61-66 NS, number of pressure steps per 
paln out 


J. 


(right justified integer).* Use -l if new 
property tables follows. 
— 2leeenumbenm Of printouts (right 


Rust ds ed) 





taeda eo weimcegers Might justified in field) 
Grominen: eG) 
Gols. 
Mo MBUs  s0 constant droplet size, =1 drop 
Peea cup 
polIwiOpem-Oc(rP) table to be supplied, =1 
X(P) is to be optimized 
13-18 MGEO, =0 circular, -=1 annular 
19-24 NDS, maximum number of S iterations 
25-30 NSO, maximum number of So iterations 
31-36 NB, maximum number of TB iterations 
37-42 NNS, first setting of step counter 
Sard 32 {format 5E12 .6) 
Co ls- 
1-12 DP, pressure step size, negative for 
decreasing, psi 
13-24 RC, mass flow ratio 
25-36 PHI, angle of annular nozzle axis, deg 
37-48 RAXO, annular nozzle axis, in. 
To-oUmrtih. toma flow rate “P final = PO + 
CNS *NP-NNS-1) *DP+DPL1 
Saeed &: (format GE12.6) 
era ssc 
1-12 H, inverse Henry's Law constant, psi 


13-24 ALAM, Lagrangian multiplier 


eo 
Ci) 





25-26 DPL, first pressure step size 
37-48 WAL, molecular weight of liquid a 
49-60 WBL, molecular weight of liquid b 
Pee (2 eo eo lener land CONnStant for component 
ge ON 

Siad io “norma: bEL2 . 6) 

eoiler 

eo Semone rland constant for component 
ory Ons 
7 Omir! drop diameter, in. 
25-36 PO, initial pressure, psia 
37-48 TGO, initial gas temperature, oR 
49-60 TLO, initial liquid tempeature, oR 
Gaewe JGO, initial gas velocity, FT/S 

Sard 6: (format 6SE12.6) 

eos). 

Pe eOeeintcials bigurd velocity, FI/S 
13-24 THOO, initial momentum thickness of 
outer wall boundary layer, in. 
25-36 THIO, initial momentum thickness of 
inner wall boundary layer, in. 
37-48 EDS, convergence criterion for §$ 
HJ=00 ESO, convergence criterion for So 
Di 2 eheeecomvergence Criterion for T, oR 
If MOP=0, the following cards are present: 


iipelOmearae’: (format 7A4., A2) cols. 





1-30 X(P) table identification (and 
alphanumeric data). 
Gaede oome(hermat 6EI 2.6 ) 
COs 
1-12 pressure, pi-2, psia 
ia edashamece. xi=2, 1N. 
25-36 pressure, pi-l So 
37-48 distance, xi-l 
49-60 pressure, pl 
61-72 distance, x1 
The last two entries are 1.0E5 and 1.0E5 right adjusted 
in their fields. The table must be monotonic increasing or 
decreasing. New property tables may be used by putting -l 
Mmomeols. 61-66 of Card 1, and following this with new 
property tables and data sets. Appendix E is a sample 


input data. 


See OUTPUT 


[meee Caso. the case identification is printed followed 
7) Se 


Meee input parameters. If MOP=0, the X(P) table forms a 
Beteenor this output. The following output then appears. 
mo xX @distaence, in. 
a »?P pressure, psia 
a R mass flow ratio 
u., vb mean free-stream velocity, ft/s 
5 oes mrew area. in. 72 
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IL 
id. 
Lae 
ol. 


14. 


nO . 
ie / . 
ns . 
ious 
a? 
Bel 


ee. 


ee 


cD 


ell 


vg 
vi 


VS 


cv 


alpha 


beta 
mg 
ro 
rog 
gonk 
wag 
wD¢ 


wg 


1b 


Posmsempe nature "Or 

liquid temperature, or 

Cas Velocity. ft/s 

igniid Velocity, tt/s 

sieeve lOClty Ve-Vl. ft/s 

Slee ract ion vs/ vb 

drop diameter, in. 

ratio of gas volume flow to liquid volume flow 
ratio of gas flow area to liquid flow area 
mass fraction of component a dissolved in 
iLatre; late 

mass fraction of component b vapor in gas 
Pace towerate., lbm/s 

imqurdes low race. Lom/s 

gas density, 1lbm/ft3 

Pacutardensiry. Lbim/ it 3 

molecular weight of component a gas 
molecular weight of component b gas 

mean molecular weight of gas 

partial pressure of component a, psia 
partial pressure of component b, psia 
latent heat of vaporization of component a, 
btu/1lbm 

latent heat of vavdorization of component b, 


muy Lem 


(en 
Oo 





seis 


Zo. 


0. 


31. 


BZ 


Bi 


34. 


oD. 


BD. 


5). 


ec . 


Bo. 


40. 


aa 


42. 


co). 


ae 


a 


oO) 


U7, 


43, 


sigma 


cgm 


em 
vigm 
vilm 
kgm 
rem 
cdm 


hm 


rb 


liquid surface tension, dyne/cm 

specific heat of gas (at midpoint of pressure 
Seo eee liom Of 

Spee ple ncaa or liquid, btu/ibm of 
WiscOcweyeOr Laqguzd, lbm/ft hr 

VEseCOsi Eyeotecas, lbm/f£t hr 

En GieieeonGgiecryilty Of gas, Dtu/hr ft of 
reynolds number of flow over drops 

drag coefficient of drops 

heat transfer coefficient of drops, 
Dae eet cD 

mass flow ratio after velocity and temperature 
equalization 

flow area after equalization, in. 2 
temperature after equalization, or 

volume flow ratio after equalization 

alpha after equalization 

beta after equalization 

gas flow rate after equalization, lbm/s 
liquid flow rate after equalization, lbm/s 
Puede secosity after equalization, 
Leimyeit. Nx 

Pas density after equalization, lbm/ft3 
liquid density after equalization, lbm/ft3 
after 


molecular weight of component a gas 


equalization 


7 





oh. 


ow . 


eu. 


a2 


a . 


oF . 


aS 


EO . 


Sy 


oo. 


oo. 


60. 


ei . 


2. 


aS . 


oo. 
BO. 


ey. 


wbgb 


web 


pab 


pbb 


ref 
yo 
wom 
tho 
delo 


delso 


redom 


cfom 


twom 
vbd 


wim 


deli 


delsi 


molecular weight of component b gas after 
equalization 

mean molecular weight of gas after equalization 
partial pressure of component a gas after 
equalization, psia 

partial pressure of component b gas after 
equalization, psia 

separator friction parameter 

separator film reynolds number 

auicediee from MOzzZle axis to outer wall, in. 
angle of outer wall relative to axis, deg 
momentum thickness of outer boundary layer, in. 
VelOegi@yeenitekiess Of Outer boundary layer, in. 
@icelaeonenestnickmess Of Outer boundary layer, 
ali oe 

reynolds number of outer Doundary layer 

Simei rietrzon ecoetficient of outer boundary 
layer 

shear stress on outer wall, psi 

mean velocity including boundary layer, ft/s 
angle of inner wall relative to axis, deg 
MOMmeheiniewentex<mess Of Inner boundary layer, in. 
TMle@ecteietnaekness Of inner boundary layer, in. 
displacement thickness of inner boundary 


Ayer. 2m. 


1 
CO 





Os . 


oa 


0). 


dig 


Ree 


1 om 


redim 


oleae 


twin 


nine 


nous 


nib 


reynolds number of inner boundary layer 

skin friction coefficient of inner boundary 
eee 

shear stress On inner wall, psi 

number of iterations required to optimize X(P) 
number of iterations required to converge cn 

Ss or so 

number of iterations required to converge on 
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IV. EXPERIMENTAL SYSTEM 


The experimental system can be grouped into three sub- 
Systems. These are: 

a) nozzle 

b) air supply system 

@) liquid injection system 
Each subsystem is described in the following sections. 


Figure 11 is an overall system schematic. 


A. NOZZLE 
The nozzle has a convergent-divergent flow passage. It 
Meml2 inches long with a variable exit’ area. The exit area 


can be varied from .45313 square inches to .84375 square 
mienes. The pivot point is located 1 inch above the throat. 
This causes the throat to vary when the exit is varied. 


Since the change in the throat is negligible, it will be con- 


Sidered to be constant. It has a throat area of .45 square 
inches. The inlet area is 1.625 square inches. The throat 
is located 4 inches from the inlet. The nozzle is constructed 


by sandwiching two 1/2 inche thick machined aluminum nozzle 
Pprorile plates between 1/2 inch plexiglas plates (Fig. 12). 
The aluminum nozzle plates are located at the end of a 30 

inch long entry section and are easily adjustable. Fi 


leas 


13 shows a close up of the aluminum section of the nozzle. 
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The inlet section provides the necessary space for pressure 
taps, liquid injection tubes, and other instrumentation and 


@ontiguration options. 


PeeeolR SYSTEM 

The air system is shown in Figure 14. Two 117 cubic 
feet tanks supplied with compressed air from an Ingersol 
Rand two-stage air compressor provides the air storage vol- 
Mm@enreguired to support the nozzle's operation. Each stor- 
age tank has a pneumatically control Norgren gate valve at 
its exit. These valves are operated by a nitrogen actuator 
and controlled by a pressure regulator. The nitrogen is 
Beeulaced to 40 psi control pressure which will open the 
Norgren valve. 

ire nozzle is supplied with air via a 3" i.d. pipe. The 
air supply to the nozzle is controlled by a solenoid actuated 
nitrogen operated 3 inch ball valve. The nitrogen is sup- 
plied via a regulator. By varying the nitrogen supply pres- 
Sure tO the ball valve, supply air pressure to the nozzle 
can be controlled. Air flow to the nozzle is measured with 
a standard ASME orifice plate. Figure 15 shows the dimen- 
emems Of the orifice. The orifice is a model D~-10512 with 


eeo.920 inch bore. 


Sees 2OULD stNJECTION SYSTEM 
@nie liquid injection system, Figure 16, is supplied by 


house water and is further pressurized by an Aurora electric 
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pump model/type 317696. The pump is rated at 50 gallons per 
min. The pump supplies pressurized water to the injection 
tube via two flowmeters (rotometers). The flowmeters are 
FEP co. precision bore flowrator tubes. One rotometer is 
rated at 1 to 12 gallons per min. and the other at .6 gallons 
per min. 

tie liquid injector is a 0.25 inch brass tube inserted 
in the 3" i.d. air supply pipe just upstream of the flange 
Bemmection to the test section. The injector tube is drilled 
with sixteen 1/16th inch diameter holes facing the test sec- 
tion entrance. The drilled holes were made as small as 
MeositbDie consistent with achieving a significant liquid 


mos tlow rate. 
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(Ta siioe noe ATTON SYSTEM 


The instrumentation system is designed for automated data 
acquisition, analysis, and display. A schematic is depicted 
in Figure 17. Where ever possible, test operations and 
memimencing are under direct computer control. Each parameter 
measured involves an appropriate transducer, excitation source, 
and calibration procedure. The major instrumentation sub- 
feyecems are: 

(A) Pressure Measuring Transducers 

M2) Nozzle Thrust Force Block 

(C) Flow Measurement Devices 


(D) Data Acquisition/Control System 


Pee RESOURE MEASURING TRANSDUCERS 

Pressure measurements are made in fourteen locations 
throughout the experimental apparatus. Eleven Micro Switch 
Meeeeeepressure transducers model PK 87633 are placed on the 
nozzle assembly to measure pressure at various axial posi- 
mbOns. Specifications for this model transducer are shown in 
Figure 18. The first pressure tap is located at one half 
miemecrom the inlet along the axis of the nozzle. The 
remainder are placed at one inch intervals toward the 


nozzle exit. These pressure taps are connected to the 


Meeeeure transducers via a 1/4" o.d. plastic tubing. The 


Oo) 





OTPBUISYOS weySAS UOTIeEJUSUINUI SUT 


~ a 


d, 


T sunsty 


WILHIANDD O/¥ 


WALSAS 
naan Like 


WI1SAS 


— WOTLIICNT ONDE 





R= eee 
Zp iviie, walle] eli < 








70 





MOUNTING DIMENSIONS (for reterence onty) 


‘ae TEM § we Tred COSTCO 


Om PrONTE 108 





140PC SPECIFICATIONS at 8.0 0.01VOC, 23°C 


ALL LISTINGS 









PRESSURE 

RANGE 
PARAMETER ee 
RGR sae oan | a ee 
a 0951 00) 168 


Eacitetion 800 200 voc 

Cutput Current 
Source 

Sink 

Supely Current Alt 

(10K ohm tosd) 















Cperaiing Temperature 


WARNING: Demage mey resuit from reversal of suppty 
and ground connections. 





meeure 18. I140PC Pressure Transducer Specifications 
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pressure transducers are mounted on the side panel of the 
Meegle. see Figure 19 for locations. Each transducer has 
an identification letter corresponding to its position. 
iiple | 1s the ID code for the transducers. 

At the nozzle inlet a Micro Switch 200PC pressure trans- 
ducer model PK 87690 is placedin the inlet line to measure 
the dual-phase mixture inlet pressure. Specifications for 
this model transducers are shown in Figure 19. Gas supply 
pressure is measured using a 200PC pressure transducer 
located at the orifice outlet. This transducer is also 
used to help measure differential pressure across the ori- 
fice. Another 200PC transducer is placed at the inlet of 
/mienorifice. The data acquisition program calculates the 
orifice differential pressure and then the air mass rate 
Meame the output from the above two transducers. 

The pressure transducers provide output voltage propor- 
tional to applied pressure. These operate from a single, 
positive supply voltage ranging from 8 to 20 VDC. The 
supply voltage was maintained at I15VDC. 

Each pressure transducer sends back a signal proportional 
to the input pressure. The signal is converted in the A/D 
converter to a digital signal which can be read by the HP9826. 
The pressure transducers were calibrated using a known pres- 
Sure source. Appendix F depicts a sample calibration pro- 


Meeiewritten for a HPY826 computer. This program reads ten 


is 








200PC SPECIFICATIONS at 8.0 t O.01VOC. 25°C (unless otherwise noted) 


Output Current Au 
oo rice 
Sink 
RY a 2) as a St ee A 
“35°C to #125° 


pereting }empature | AU 
| (-65° F to +257° F) 





°F $.0. 2 he argedracc Gflerence Detween end pants (null and full prespure oulpum) 


Masimum pressure wi(houd Carrage 
Wrinoutl housing envelope rupture. 


ELECTRICAL AND PRESSURE CONNECTIONS 





meeure 19. 200FC Pressure Transducer Sp 
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faeces! Loe code for Instrumentation 


Meeter Designation Description 
A HOO from 2tnlet 
B ao! eon Inlet 
G eo eeeoOm. 11 Let 
D oy! Eromoinlet 
iB Beso enmom- 2 Nike: 
i oo 'Memom inlet 
G To  Serom inlet 
A fn) wecom anler 
if Po Senom Ln vet 
J foo oem ime c 
K Oo semen Ther 
L inlet nozzle pressure 
M force-block signal 
N force-block excitation voltage 
0 orifice exit pressure 
PP orifice inlet pressure 
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Pressure = A + B(volts) + C(vol ts)? + D(volts)? 


Pressure 


transducer A B C 0 
A “6.049 3.926 -.17336 «> .01121 
B -7.059 4.234 -.22921 01415 
C -.6491 3.934 -.16025 00945 
D -7.056 4.692 -. 34118 .01898 
is _27.009 4.152 -.20717 .01237 
F -6.389 3.933 -.16055 .00957 
G -6.589 3.956 -. 16427 .00961 
H -§.588 3.908 -. 15806 .00942 
I -6.595 3.949 -.16270 .00957 
J =Oecod dade? -. }6509 .00998 
K -6.203 3.846 -.14697 .00880 

eu -20.207. 11.157 -.08591 00492 
0 -20.668 11.068 -.07562 .01118 
P “18.827 11.928 -.17715 01101 

igure 20 bese Ob FOolynomial Coerficients 


for Pressure Transducer 


lis 





values of pressure for each 140PC transducer and gives the 
mean and standard deviation for each from 0 to 45 psi in 5 
Bei increments. Appendix G depicts sample output of the 
program for a given pressure. Appendix H illustrates the 
Meeeram used to calibrate the 200PC pressure transducers. 
The program works as above, except values are taken from 0 
to 60 psi at 5 psi increments. 

Data obtained during calibration source pressure lis 
plotted for each pressure transducer (Appendix I). Each 
plot was curve fitted with a third order polynomial. Figure 
20 shows coefficients of the polynomial. These polynomials 
are used in the data acquisition/control program to convert 


transducer readings to pressure readings. 


Bee NOZZLE THRUST FORCE-BLOCK 

the thrust produced by the nozzle was employed to deduce 
the nozzle exit velocities. The thrust was measured by in- 
strumenting a target plate in the exit flow field. Appro- 
priate screens were installed to prevent liquid "bounce 
Dack." The jet momentum force on the target 1s acquired by 
a Dalance beam system shown in Figure 21. A Kistler-Morse 
force block provides an analog signal proportional to the 
nozzle jet momentum. 

The calibration of the force-block was completed by 


placing known weights on the force-block side of the balance 
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beam, (Figure 21), and recording the voltage produced. 
Appendix J shows results of this calibration. These results 
were curve fitted with a third order polynomial. The follow- 
mieis the volt-—to-force conversion polynomial: 
FB = (1/453.6)*(-234 .3964+48715 .28*VOLTS~11261 .54*VOLTS+33899 . 3*VOLTS) 
mie FB“LB/LN 


FN POR@e ALT NOZZLEVExXIT CLBF) 


mem- FORCE AT FORCE BLOCK (LBS) 
mee= LENGTH FROM PIVOT TO FORCE BLOCKCIN.) 


mie- LENGTH FROM PIVOT TO DIRECTIONAL OBSTRUCTIONCIN.) 


eee FLOW MEASUREMENT 

1. Air mass flow measurement for the system achieved 
by measuring the inlet and outlet pressure on a D-10512 ori- 
Meee with a 0.92 inch bore. For details on orifice see 
Meeeron Tit. Air flow is calculated in the data acqui- 
SeeedoOn program. The program uses the following equation 


Obtained from References 4 and 5. 


= ] i Tan ee oes ey : 
M . i A, 26 p,¢P1 P2 where 


ea 


B= CE 
2 Peg yl! 2 


m= d/D 

Bee 0.60 

A,= AREA OF ORIFICE = mD* FT? 
a 





cm eee / SEC 2 

oie BENSe iy CLBM/ PLS 

Ree ReoSoURE AT INLET IN LBF/FT2 
Pemeee ee cooURE, Al LUTLET IN LBF/FT2 
en ATR MASS FLOW RATE LBM/SEC 


ime discharge coefficient C is the factor that accounts 
for losses through the orifice. Since the values of C varies 


T with B = .3, 


maem .02 to .60 for Ra number from 1" ee ILO 
Mew2ll be considered constant. 

2. Water flow measurement was made using two rotometers. 
Calibration of the rotometers were made by measuring the 
time for a given weight of fluid flow. The mass flow rate 
was calculated and plotted versus the rotometer reading. 
Appendix K is the plot of the results. The plot was curve 
meered With a third order polynomial. The following is 


that polynomial: 


Z 2 


Meo = —.0063268+ .002097278"*RR-.00000658*RR°+1.1X10*RR 


where RR = rotometer reading 


See tA ACQUISITION/CONTROL SYSTEM 

The heart of the data acquisition system is a HP9826 
Smet computer. The HP?825 communicates via a Hewlett 
Seer d 3497A data acquisition/control system. This system 
gathers data from the pressure transducers and nozzle thrust 
Memee—block. It converts the analog signal to digital data, 
and stores the data in memory. Figure 22 shows a pressure 


Peemsducer to A/D converter channel connection, and Figure 23 
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the force-block to A/D converter channel connection. Figure 
Meedepicts the current connection to the data acquisition/ 
Sonerol system bus from the sensors. 

The data acquisition and analysis program is a fully 
automated system which integrates the use of the HP9826, 
HP3497A, and all sensor devices. The program can be divided 
into three sections: 

fee Data Acquisition - Data is acquired from 15 differ- 
ent sensor devices. All system parameters are ini- 
tialized and configuration types are inputted into 
the computer memory for future use. Ten readings 
are made on each sensor and a mean value for each 
is calculated and stored in memory. 

2. Data Conversion - Input data conversion is completed 
by using calibration information obtained on each 
pressure transducer and force-block. The program 
takes the stored mean values from memory, and 
Poem verts the 2nput voltage to a pressure or force 
Heameowene third order calibration polynomial for 
each sensor. 

3. Data Analysis - This is the heart of the program 
in that all performance parameters are calculated 
and printed. It uses all of the above information 
to calculate mass flow rate of water, mass flow 


Pageowordir, tOLat mass flow rate, thrust, mixture 





ratio, exit velocity and all pressure informa- 
Zion taenrOllowing equations are also used: 
Exit Velocity = (Thrust/Total Mass Flow Rate)* g (ft/sec) 


/m 


m : 
water air 


Peaexture Ratio 


ee, (ibm 


Total Mass Flow Rate = Mater air ‘sec 


Appendix L is the data acquisition program and 


Eample output. 
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eee SEN OAL RESULTS 


Six test runs were made using the equipment described in 
the previous sections. It was found that the inlet pressure 
to the nozzle remained nearly constant for each setting. 

The testing was conducted by maintaining a constant air 
control valve setting and varying the mass flow rate of 
water. The regulator outlet pressure of the nitrogen bottle 
which controls the air valve was set at 30 psi and 15 psi. 
These settings correspond to nozzle inlet pressure of 35 + 
Meeepsi and 29 + 1.0 psi respectively. For each nozzle in- 
let pressure, experimental data was obtained at varying 
mass ratios of water-to-air in the range of 2 to 13. These 
experiments were conducted with three discharge area ratios 
ot Aoyt): Those ratios are: Ceo, wee O0 .- and 
1.9259. The results of the six experimental tests are dis- 
@eayed in Tables II through VII. 

Inlet air was limited to a maximum value of 40 psi. 

The 140PC pressure transducer's maximum output pressure is 
Sempsi, (See information on 140PC transducer Figure 21). 
The experimental system has the ability to reach higher 


Meeooures. The present configuration may be operated to 


mo) psig. 


re 
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Vit=. @DESCUSSION 


The results of this investigation are in two parts; the 
experimental test results and the computer outputs based on 
initial conditions similar to those of the experiments. 

The key variables are: the nozzle overall area ratio 
Ags Brquid/air mass flow ratio, nozzle exit velocity, nozzle 
supply pressure, nozZle thrust, and the nozzle axiai pres- 
Sure profiles. Figures 25 through 34 present the exit vel- 
ocity vs. mass ratio results for the experiment and computer 
output. Figure 29-34 illustrates the comparison of the ex- 
perimental results and the corresponding computer output. 
Meeure 35 illustrates the variation of nozzle thrust with 
mass ratio for different area ratios and inlet pressures. 
Figures 36 and 37 present the axial pressure profiles. 

In all cases the exit plane velocity decreases as the 
liquid/air mass ratio is increased. In the low mass ratio 
range (1.e., less than ~ 5) the velocity decrease is very 
pronounced. Past a mass ratio value of ~ 10 the velocity 
of the mixture becomes relatively insensitive to the mixture 
ratio. This is as would be expected. For the same inlet 
conditions there is a fixed amount of energy available for 
Senversion to its kinetic form. In the nozzle process the 


energy is conserved and thus an increasing mass ratio is 


a2 





manifested in a decreasing exit plane velocity. In all cases 
the agreement between the experimental tests and the analytic 
Seedictions are with in ~ 10%. 

There also was apparent and a consistent trend with re- 
spect to the nozzle area ratio and the exit velocity. As 
the nozzle exit area was increased the exit velocity of the 
Maxeture decreased. The effect seems to be more pronounced 
at the lower mass ratios. This trend was confirmed by ob- 
servations of the actual flow field in the nozzle. If the 
exit area was increased significantly (beyond the max area 
used in those tests) there was evident an abrupt and severe 
separation within the diverging portion of the nozzle 
passage. This was accompanied by a drastic decrease in 
exit velocity as evidenced by the output from the thrust 
target. It appears that the diverging portion of the 
nozzle, at a certain point, starts to behave as a subsonic 
diffuser and hence experiences an adverse pressure gradient. 
This reasoning is in part confirmed by the pressure profiles 
meencictied in Figures 36 and 37. 

The relationship between measured thrust and mass ratio 
(Figure 35) indicates a slightly increasing trend. This may 
De explained by considering the following. Each set of data 
points (Cat constant inlet air pressure and constant nozzle 
€xit area) is developed by varying the mass ratio. This, 
in turn, is achieved by increasing the liquid rate by in- 


Sreasing the liquid supply pressure. The net result is an 
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increase in the liquid inlet velocity or an increase in the 
inlet energy level. Thus a particular data set is not truly 
at a constant inlet energy level but is increasing. 

Perhaps the major discrepancy between the experimental 
test results and the computer analytic model is evidenced 
in the axial pressure profiles of Figures 36 and 37. [It is 
Silearly evident that at a certain point the nozzle passage 
reverts from a nozzle to a diffuser. This transition point 
Beeurs slightly downstream of the throat and the pressure 
starts to increase. 

Unfortunately the computer analytic model requires a 
pressure profile as an input. Furthermore the pressure pro- 
Bebe must be continuously decreasing. Thus the profile as 
obtained from the experiment are not directly useable. The 
Meedation is examined with the aid of Figure 38. Curve A ls 
a typical axial pressure distribution as obtained from the 
experiment. Curve B is a pressure distribution used by 
Elliott in the application of his computer program. Pressure 
profiles C and D were arbitrarily defined and the exit plane 
velocity for each was calculated. Velocity variation in the 
Pange of 1% is evident. it appears that the final exit 
velocity is relatively insensitive to the actual pressure 
Deorile in the nozzle. 

In light of the preceeding difficulty of matching the 


experimental pressure profile to the computer model and the 
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apparent flexibility of the type of pressure profile it was 
decided to employ a form of profile D of Figure 38. Thus 
instead of using the actual experimental pressure data as 
an input in the computer model, a profile resembling curve 
D was developed for each inlet pressure case. 

It 1S apparent that within certain rather wide limits 
the Elliott computer model yields results that correspond 
Within 10% to results obtained from the experiment. The 
general trends have been confirmed and their behavior has 
revealed nothing unexpected. The conclusion of Elliott 
[Ref. 1] has been largely substantiated. "It is very 


difficult to design a bad or a good dual-phase nozzle." 
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WOOD BELOIT DEIONS IEG) 


A set of experiments were conducted with an air-water 
mixture flowing at a mass ratio from 2 to 13. The experi- 
ments correlated well with a dual-phase two-component 
computer program. It appears that the program will permit 
prediction of the exit velocities to an accuracy of 10% 
for two-phase nozzles. At small nozzle area ratios, accu- 
racies between theoretical and experimental are better; 
approximately 5%. Ail predicted velocities are higher 
than measured experimental values due to estimation of the 
drop size and the initial kinetic energy of the liquid at 
mayection point. 

It would be desirable to develop a drop size subroutine 
to better estimate varying drop size. The experimental 
system can be improved by velocity measurement devices at 
the inlet. Input of the gas and liquid velocities are 
necessary in the dual-phase two-component program and thus 
these measurement devices are vital for better accuracy. 

Due to the insensitivity of the pressure distribution, 
it appears that any reasonable approximation to the pressure 
profile can be employed in the two-phase two-component flow 
Program. For a given nozzle exit area and inlet pressure, 
a pressure distribution curve can be obtained through the 
experimental system. This distribution can therefore be 
used to obtain good results. 
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SAMPLE P(X) PROGRAM 
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SAMPLE HEAT CAPACITY PROGRAM 
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PPPENDIX IT: CALIBRATION PLOT OF PRESSURE TRANSDUCERS 
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APPENDIX L: DATA ACQUISITION AND ANALYSIS PROGRAM 


S PRINT “FUT AGTOMETER READING” 
4 INPUT Rr 

5 PRINT Rr 

Fi DIM A(14),B(14) ,C(14) .Dist(14) ,D(14) 
8 DIM Press(14) 

10 DIM X(18,100),.Sum(18) ,Mean(18) ,Sum2(18) 
30 N2=1Q 

40 FOR IJ=1 TO N 

S50 SUTrUT e703: "ALIVii" 
60 SHER dee X (1, 1} 

70 QUTPUT 709;"AI2VT1" 
80 SNIER #70, Xt...) 

30 OUTPUT 70S9:"AISVT1" 
100 ENTER 7093X(3.1) 

110 OUTPUT 709; "“AL4VT1" 
12 ENTER 709:X(4,1) 

130 OUTPUT 709; "AISVT1" 
140 ENTER. 2093 %05,1) 

150 OUTPUT 709; "AI6VT1" 
160 Sen 7 ue. CG. > 

170 OUTPUT 709;"°AI/7VT1" 
180 ENTER 7093X(7,1) 

130 OUTPUT 709; "AI8VT1" 
200 ENTER 709:X(8,1) 

210 OUTPUT 709;"AI9VT1° 
220 ENTER 7092 4(9..1) 

230 OUTPUT 709: "AI10VT1" 
240 BIER 703.401 05 1) 

250 OUTPUT 709; "AI11VT1" 
260 ENE  /US34 C11, 1) 

261 SU eu 709-sAli2ZVii* 
262 Bier 70s Gl 12.0) 

263 GUTPUT 709;"“AI13VT1)" 
264 ENTER 709:X(13,1]) 

265 GUTPUT 70S: "Al14VT1)" 
266 ENTER 709:3X(14,1) 

267 QUOTES 70S: AlitSvt1)” 
268 ENTER 9709-4115, 1) 

2/0 NEXT I 

440 FOR J=1 TO 16 

450 Sum(J) =0 

460 FOR I=! TON 

470 Sum(J)=Sum(J)+X(J. 70) 
480 NEXT IT _- 

4930 Mean(J)=Sum(J)/N 

510 NEXT J 
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DUAL-PHASE TWO-COMPONENT COMPUTER PROGRAM 
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CRITE (6,9352) 

EP (kTBgT} 403 403,400 

READ (18) T3PRAYI 

@RITE (6,9491) (TORAYI(T),I=1, 100) 

KTSCTSKTECT~ 1 

GO TO 399 

RITZ (60,9491) (TBRAY(I) ,L=1,KTB) 

EAR TB SAVE ARRAY. 

XT B=1 

KTSCT=0 

DO 305 r=1,100 

TBRAY (I) ="0.0 

CONTINU 

LL=25 

GO To 500 

ORE NUMBER OF ITERATIONS REQUIRED TO CUNVERGE ON TB. 

2eeeeee Ce O88 OE & STATES NT NG. 140 SOR eSSCHHE REET AHK EASE KRESSEHKESe HHH EEE 
Seeeee ee AeReae ee CaSe eae SSeAeee es CSHaeeee SHER ARBSEAEABESERERESEAEBEHEEKLEASEEASE EE 
NTB=NND 

NNB= 


MPUTE DENSITY OF GAS MIXTORE. 

RHOSB= (WGB*P) / (10.732*TB) 

TERPOLATE POR DENSITIES OP LIQUIDS A AND 3. 
N= 10 


NNe ha} 
CALL INTEP (ROALB,N, TB, P) 
Gue27 
Mais) 282s 319,360 
¥= 
NN=11 
ee INTR? (2OS8LB,N,T3, P) 
PEN) 3 19,319,363 
GCAPUTE OENSITY OF LIQUID SIXTURE. 
lS Re A hue Bal ae - O- ALP OLS) AVR OBLB) )) 
See *ENGP)*® (1 HOG B+ 2B/ RHO La) ) 
RVE=SEOQOLE/S Op eegar ot 
Q= (QHOLB*ABAR* VB) /(144.0*E8LB) 
NTERPOLATE INTO VISCOSITY TABLES POR LIQUIO A AND 8. 
N21 
NNW=1% 
SAbGel YT SPAY LALO, dy lO, ©) 
Lis2d 
tres 19.319, 366 
N=13 
NN=15 
CALL INTEP(VIHLB,N, TS, P} 
IF (s) 369,309,370 
LL=28 
GO TO 319 
re Re ee OAL EE ea taty 
Torey Veo ce boca) /iviLa® serege end 
ER SGUCS ENT VALUGCS OF A AND ?P NTD A475 27 TADLA 
NSaess929 
RETURN 
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